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Abstract 

Recent experiments have demonstrated that the glass transition temperature of thin polymer 
films can be shifted as compared to the same polymer in the bulk, the amplitude and the sign 
of this effect depending on the interaction between the polymer and the substrate. A model has 
been proposed recently for explaining these effects in two limiting cases: suspended films and 
strongly adsorbed films. We extend here this model for describing the cross-over between these 
two situations. We show here how, by adjusting the strength of the adsorption, one can control 
the glass transition temperature of thin polymer films. In particular, we show that the shift of 
glass transition temperature, refered to that of a suspended film, varies like e 1 / 72 where e is the 
adsorption energy per monomer and 72 is the critical exponent for the mass of aggregates in the 
2D percolation problem. We show also that the interaction leads to an increase of T g for e as 
small as a few 0.01T where T is the thermal energy, and that this increase saturates at the value 
corresponding to strongly interacting films for adsorption energies of order 0.5T. 

PACS numbers: 64.70.Pf Glass transitions, 61.41.+e Polymers, elastomers and plastics and 68.15.+e Liquid 
thin films 
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I. INTRODUCTION 



When cooling a glass forming liquid, instead of freezing at some well defined temperature, 
one observes a huge increase of the viscosity which takes place continuously. Such glass 
formers can be either simple liquids or polymer liquids, and many features are similar in 
both regarding the glass transition. One defines the glass transition temperature T g as the 
temperature at which the dominant relaxation time on the molecular scale (or monomeric 
scale in the case of polymers) reaches about ~ 100s, which corresponds typically to a 
viscosity of 10 12 Pa.s in the case of simple liquids. Typically for such glass forming liquids, 
the viscosity increases by 12 orders of magnitude over a change of temperature of about 100i^ 
down to T g . For reviews on the glass transition see e.g. The underlying mechanisms 

involved in this dramatic increase are still poorly understood. One of the issue has long 
been to know whether the relaxation mechanisms take place on an increasing length scale 
when lowering the temperature Q. This issue .ead to a huge interest in studying the g.ass 
transition in thin liquid films over the past ten years [4J]. For practical reasons (film stability) 
these studies have been performed on thin polymer films. Indeed it has been demonstrated 
that the glass transition behaviour in thin polymer films is very different from that of the 
bulk. When considering a film deposited on a substrate with which the interactions are 
supposed to be very weak, Jones et al [5( have shown that the glass transition temperature 
can be reduced by as much as 25K for 100A thick films. Numerous studies have confirmed 
this effect since then (see e.g. 00630). In the ease of polymers with radius of gyrations 
larger than the thickness of 7he film, Forrest and co-workers Q 3 found even larger T g 
reductions, by as much as 70—80 K for 10 nm thick films. In this manuscript we will consider 
only the case of shorter chains, with radius of gyration smaller than the thickness of the film. 



On the contrary many other experiments performed on str ongly a dsorbed fi 
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ms have shown 
, sometimes by 



that the glass transition temperature increases in this case 
as much as 60K for films 80A thick. More generally, it has been demonstrated by various 
authors that both the sign and the amplitude of the shift depends on the interaction with 
the substrate 

0, H, 0. These experiments were performed by measuring the 



thermal expansivity of the films by ellipsometry 
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or Brillouin 
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light scattering [6j], by probing molecular dynamics using dielectric spectroscopy 



by calorimetric measurements 
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Note that these T g shifts are very important effects if one considers molecular relaxation 
times, or viscoelastic properties. For instance -according to the usual WLF law- a shift of 
T g of a few tens of Kelvins results in an increase (or decrease) of relaxation times by several 
orders of magnitude at a given temperature T. The same polymer film is therefore expected 
to have fundamentally different viscoelastic properties depending on the substrate on which 
it is deposited. Note however that most measurements of T g shifts have been performed either 
by dilatometry (ellipsometry) or by dielectric measurements and are somehow indirect. More 



direct measurements by Atomic Force Microscopy have been reported e.g. in 



22, m 



but brought conflicting results. Indeed, some studies have confirmed previous results [2J 



231 ] while another 2l( failed to reveal corresponding changes in the viscoelastic properties. 



However, even though experimental studies are still required to get a complete understanding 
of thin film behaviour, it is generally considered now that huge dynamical changes take place 
in thin films as compared to the bulk. 

Another feature of glass transition which has been emphasized recently is the stro ngly 
toerogeneous nature of the dynamics close to T s QQQBHQHBHQQfl 
1361 ] . This heterogeneous nature has been demonstrated by many experimental studies, using 



;echniques such as NMR 



31 



Q 33,34,35, 3(J, 



24, 



25l l26j]. fluorescence recovery after photo-bleaching (FRAP) 



translational and rotational probe diffusion 



burning [27] or solvation dynamics 



dielectric hole- 



]. These studies have demonstrated the coexistence 
of domains with relaxation time distributions spread over more than 3 to 4 decades at 
temperatures typically 20 K above T g . The characteristic size £ of these domains, measured 
by NMR [2^1, is typically 2 to 4 nm in the case of van der Waals liquids and as small 
as 1 nm in glycerol j^J. Using a deep photo-bleaching technique, Ediger and co-workers 
have measured the life-time of these heterogeneities, which they found comparable and even 
larger than the dominant relaxation time t v at temperatures below T g , while Spiess 
and coworkers j^l found both times to be comparable. For reviews on these issues, see e.g. 

y y Bsi- 

One of us recently proposed a model for the glass transition in the bulk, which allowed 
both to interpret thin film experiments, and to prop ose an explanation for the origin of the 
heterogeneous nature of the dynamics in the bulk |40ll41|. in the case of van der Waals liquids. 
We proposed that the glass transition in the bulk is controlled by the percolation of slow 
domains corresponding to upwards density fluctuations on the scale of a few nanometers. 
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Note that the possible role of density fluctuations had first been considered in references 



42, 



43( and that the model proposed by Long et al can be considered as an extension of 



that proposed by Ediger. In the case of a suspended film, the glass transition occurs when 
the slow domains percolate in the direction parallel to the film, which requires a larger 
fraction of slow domains and results in a decrease of the glass transition temperature. In 
the case of films deposited on a substrate with which the interactions are strong, the glass 
transition occurs when both interfaces are connected with continuous paths of slow domains. 
This requires that the correlation length of the 3D percolation problem is comparable to 
the thickness of the film, which requires a smaller fraction of slow subunits as compared to 
the bulk glass transition. In both cases -suspended films and strongly interacting films- we 
predict a shift of the glass transition of the form 

^ = <) 1M (1) 

where (3 is a number of order unity. The sign of (3 is positive in the case of strongly 
interacting films (SI films) and negative in the case of suspended or weakly interacting 
films (WI films). z/ 3 is the critical exponent for the correlation length in 3D percolation. 
Note that the mechanism at work in the model proposed by Long and co-workers applies 
only for simple liquids, or short polymers, in that it does not predict molecular weight 
dependence of the T g shift such as those observed by Forrest and co-workers 0,0). Thereby 
we restrict our discussions below to the case of low molecular weight polymers, i.e. chains 
with radius of gyration smaller than the thickness of the film. Note that de Gennes 0] has 
proposed a mechanism at work in the case of large chains, which is not incompatible with 
the mechanisms described in j^flj. 

The aim of this paper is to describe the cross-over between the two extreme cases, and 
thereby to give a more precise meaning to the expressions "weakly interacting films" or 
"strongly interacting films". The issue is: what is the strength of the polymer-substrate 
interaction required for observing an increase of the glass transition? Are typical adsorption 
interactions sufficient or is chemical grafting required? The paper is organized as follows. 
First (section II) we summarize briefly the model for the glass transition in van der Waals 
liquids proposed in j^J. Then in section III we re-discuss the glass transition in thin 
suspended films or in strongly interacting films. In section IV we discuss the cross-over 
between the WI case and the SI case precisely, and provide an expression of the energy scale 
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of the interaction at the cross-over. We discuss our predictions (section V) and compare 
them to some recent experimental data (section VI). Note that in a joint paper [45], we 
study in more details percolation properties in thin films in various situations corresponding 
to the physical cases studied here, and we perform Monte Carlo numerical simulations to 
illustrate and support our theoretical considerations. 



II. MODEL FOR THE GLASS TRANSITION IN VAN DER WAALS LIQUID 

As mentioned in the introduction, it has been demonstrated over the pas t ten years that 



;he dynamics is strongly heterogeneous close to the glass transition 
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38, 



39|. One of us 



4l| proposed that the slow subunits correspond to upwards density fluctuations, and 
fast subunits to downwards density fluctuations. We introduced a density p c above which a 
microscopic domain of volume Vq (to be defined below) is in a very high viscosity state, i.e. 
has a relaxation time of order r g , where e.g. r 9 ~ 100s is the time scale chosen for defining the 
glass transition. We proposed that the scale vq is determined by two competing relaxation 
processes: individual monomer jumps and the relaxation of the density of dense subunits 
in the more mobile neighbouring subunits by a diffusion process. The competition of these 
two processes results in the building up of the slowest relaxation process at an intermediate 
scale that we calculated to correspond to a number N c of about a few hundreds of monomers 
(4]], which corresponds to a characteristic length £ of 2 to 4 nanometers and to a volume 
Vq ~ N c a 3 of order a few 10nm 3 typically (a ph 5 A is a typical monomer length), in agreement 
with experimental results. The probability distribution of the density of these small domains 
of volume Vq is 

P(T, P )~e xp C M 2r 2 P)2V ° K ) P) 

Al Peq 

where p e q{T) is the mean equilibrium density at temperature T and K is the bulk modulus. 
In our model, the glass transition is defined as the temperature at which subunits of volume 
Vq and density larger than p c percolate. According to our model, the heterogeneous nature 
of the dynamics results from the strongly non-linear behavior of the dynamics as a function 
of p close to T g . More generally, at any temperature T, the viscosity is dominated by the 
characteristic time scale corresponding to the density p c (T), such that subunits of density 
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equal or larger than p c (T) percolate, that is 

P(T,p)dp = pl D (3) 



Pc(T) 



which reads 



Pc(T) p eg (T) = JL)i/* F ( p 3D) ( 4 ) 

Peg KV 

where F(x) is the reciprocal function of the error function Erf(x) defined by 

1 f + °° -u 2 
Er f( x ) = , n m/9 / exp(— -)ciu (5) 



(27T) 1 / 

Note that F(x) is positive for x < 0.5 and negative for x > 0.5. pl° is a 3D perco- 
lation threshold. Thereby, the dominant relaxation time in a mechanical experiment is 
t = T(p r (T)). It coexists with much faster relaxation times which are not relevant for the 



viscosity [40j,|4l|. Knowing N c , one can calculate p c — p eg (T g ), using equation (HJ). If one 
considers for instance bond percolation on a FCC cubic lattice, one has p^ D ~ 0.119. With 
the typical values N c ps 100 and K ~ 10 9 Pa, one obtains that (p c — p eq )/Peq is of order a 
few percent. 



III. GLASS TRANSITION IN THIN SUSPENDED FILMS AND IN FILMS WITH 
STRONG INTERACTIONS WITH THEIR SUBSTRATE 

A. Glass transition in thin suspended films 

Consider first a film of small thickness h of order a few hundred Angstroms, either freely 
suspended or deposited on a substrate with which interactions are weak. In the same way 
as for the glass transition in the bulk, for the thin polymer film to be in the glassy state, 
we assume that the viscous domains of volume v o and density larger than p c percolate. The 
percolation that is required is in the direction of the plane parallel to the film so as to build 
macroscopic clusters of slow dynamics. On the scale £ = aN^ 3 , the film has a finite number 
n = h/l; of layers in the direction normal to the film. For £ w 2nm (which amounts to taking 
N c w 100), n is of order a few tens. Then, as far as percolation is concerned, the situation 
is intermediate between a 3D and a 2D problem: we call it a quasi 2D problem. One can 
show j3| that the percolation threshold p^(n) in a film with a finite number n of layers is 
larger than the 3D percolation threshold and is given asymptotically by 

p+(n) = p 3 c D + p n - 1/v3 (6) 
6 



where \x is a number of order unity (~ 0.5 for a cubic lattice), and z/3 « 0.88 is the critical 
exponent of the correlation length in 3D percolation. Thus, the cross-over between 2D and 
3D is described by a power law and is wide. For the film to be in the glassy state, one needs 
a fraction p^(n) of slow sub units, which is larger than that required in the bulk, p\P . By 
calculations similar as those used to calculate p c , it has been shown in j3] that the relation 
between the temperature T and the fraction p of subunits of volume vq with density larger 
than p c is given by 



where we have used t 
expansion coefficient) 
can be linearised to yield 



l ^ = (F(p)-F(pl D ))-^ (7) 

ie relations (K / '{T 'p eq )) l ^ 2 XT ~ 1 (xt is the dimensionless thermal 
Since we consider only small differences for p— p z c D , this equation 



3D T-T s Nl p T-T, 

with F'(p) = dF/dp and 7 = —N^ 2 /F'(pl D ) > 0. Writing equation (jSJ) with p = p£(n), 
and using equation ©, one obtains the glass transition temperature T~(n) of the film: 

= y fef (n) " pT) = ^n-^ = ^{\)~^ (9) 

Note that we are using the notation p^(n) here because p^{n) is larger than p z c D and AT~(n) 
because T g is shifted downwards in this case (AT~(n) < 0), as it is clear in equation Q. 
From this result, one deduces the T g shift: 

Tg ^(1/2-1/(3,3))^^ P V 1 J 

which we can write also 

T^(h)=T g (l-pfy 1 ''*) (11) 

with (3 = — pF' (pl D ) / Nc 1 ^ 2-1 ^ 3 ^ and a is the monomer length (typically 5A). With N c ~ 
100, v « 0.88, n « 0.5 Q, Pe D = °- 119 H and ^'(0.119) « -4, one obtains that (3 « 1.5. 
In the case of polystyrene, for which T 9 w 375K, this leads to a decrease of T g of about 20K 
for a film 10 nm thick, which is in agreement with typical experimental results 0,0|. One 
of the most important conclusion here is that the T g shift that we predict is long ranged. 
According to our model, the T g shifts measured e.g. in ellipsometric experiments such as in 
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references p, |8J do not result from an average effect of strongly perturbed thin layers at the 
interfaces of the film: the "bulk" of the film itself has a lower T g , as given by equation ([lip. 

Fig.l illustrates our model for the glass transition in thin suspended films. One can see 
from above a spherical body, with a radius larger than the thickness h of the film. Both 
caps of the particle emerge from the film on each side. At a temperature higher than 
T g + AT~{h) the slow aggregates of relaxation time r g may connect both surfaces, but they 
do not percolate in the direction parallel to the film. They build only finite clusters, which 
move freely in the film. The viscosity of the film is measured by pulling the spherical body. 
It is smaller here than the viscosity in the bulk and it is dominated by the slowest subunits 
which percolate in the direction parallel to the film. The measured viscosity r}fu m of the film 
may be calculated from the T g shift that we calculated before : 

r] fil m(T, h) = r] bulk {T - AT g -{h)) (12) 

where AT~{h) is given by equation (|TUj) and 775^*. (T) is given by the corresponding WLF 
(or VFT) law of the considered liquid (remember that T — AT~(K) is higher than T since 
AT~(h) is negative). The viscosity which dominates the motion of the body is that in the 
middle of the film, where it is the highest. Note that in a two layer model, in which the 
measured T g shift would result from an average between a fluid thin layer and an unperturbed 
"bulk" of the film, the viscosity of the film would be given by the viscosity of the bulk. 
Thereby, an experiment analogous to that proposed here (cf Fig.l), or an analogous one, 
is discriminating between our interpretation and that of the two layer model j?l Isf • At the 
temperature T g + AT~(h) (Fig. l.b.) the slow subunits percolate in the direction parallel 
to the film and the film is at its glass transition: in order for the body to move, the slow 
subunits have to deform and reorganise. The viscosity of the film measured by pulling the 
body is thereby the same as that of the bulk at T g . 



B. Glass transition in films with strong interactions with their substrates 

Let us consider now a thin polymer film with strong interactions with the substrate, at 
a temperature T > T g (bulk). The correlation length of slow subunits, i.e. the size of slow 
aggregates is ( ~ £|P — Pc D \~ V3 -where p is the fraction of slow domains at temperature T. If 
p is sufficiently close to pl°, i.e. if the temperature is not too far above T g , the size of slow 



S 



aggregates is comparable to the thickness h of the film. This means that there is a finite 
fraction of monomers from the free interface which are connected to the rigid substrate 
by continuous path of slow subunits. Then, the film is glassy. Indeed, consider a shear 
experiment between the two interfaces: there is a no-slip boundary condition on the upper 
interface, while the monomers are attached to the other: one needs therefore to deform slow 
aggregates which thereby control the viscosity of the film. The condition for being glassy 
reads then: 

c-eb-pf r s ~h (13) 

Thereby a fraction of slow subunits p~ (n) smaller than pj? D is required at the glass transition. 
It is given by 

pf -pjW^n- 1 ^ (14) 
which results in an increase of T g given by 

AW „„ a 

which we can write also: 



lli = p>qy/>* (is) 



. a , 



W=T,(l + / 3'(p lM ) (16) 

where f3' m —F'(p i c D )/Nc 1 ^ 2 ~ 1 ^' iu ' i ^ is analogous to (3, with a different pre-factor however. To 
illustrate our model, one may also adopt an equivalent point of view, by considering a semi- 
infinite polymer melt in contact with a strongly interacting substrate, at a temperature 
T > T g where T g is the glass transition temperature in the bulk. Monomers directly in 
contact with the substrate are immobilized. Then, one may ask: is the melt in a glassy 
state in the vicinity of the substrate? If yes, what is the range of this effect? We argue 
that at a distance z smaller or comparable to the 3-D percolation correlation length ( ~ 
i\p~p1 D \ 1/3 -where p is the fraction of slow domains at the temperature T-, a finite fraction 
of monomers are connected to the substrate by continuous paths of slow dynamics domains. 
Then, these monomers will control the dynamics at the distance z, as a consequence of the 
strong adsorption. Indeed, without adsorption, the continuous path may move freely within 
the melt since T > T g . Then, the polymeric liquid is glassy within a distance d to the 
substrate given by 

d = b'(-^-r (17) 

1 1 9 
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with b' = a(3' U3 is of order one monomer length, i.e. 5A Note that at about T = T g + 10fC, 
the distance d can be a few tens of nanometers. We have just emphasized that the 
effects we predict are long ranged. Indeed our model predicts a shift of T g given by 
a power law as a function of the distance to the interface. In particular, we argue 
against the interpretation proposed by some authors (see e.g. according to which the 
measured T g shifts result from an average effect, between strongly perturbed thin interfacial 
layers and an essentially unperturbed "bulk" o 
one considers viscoelastic properties of thin films 



thin films. This point is essential when 
or of filled materials 



21 



22 



47 



48j,|49|. 



The case of a SI film is schematised on Fig. 2. Let us consider here a body, such as the tip 
of an AFM, penetrating in the film, within a distance z to the substrate. We suppose that 
the interactions between the tip and the polymeric liquid are weak, whereas they are strong 
between the liquid and the substrate. There is indeed a no slip boundary condition on the 
surface of the tip, whereas the monomers in the aggregates are attached to the substrate. 
At a temperature T larger than T g , and at a distance z larger than the size ( of the slow 
aggregates, the viscosity controlling the motion of the AFM tip is smaller than that at T g . 
Indeed, at this distance, the slow aggregates move freely in the film and can move away 
rigidly to open a path for the tip. When the tip gets closer to the substrate, the viscosity 
increases (Fig. 2.b). Indeed, at distances smaller than (, the AFM tip has to deform the 
slow aggregates in order to move parallel to the film. More precisely, the viscosity measured 
by the probe is then: 

r)fiim(T, z) = Vbulk {T - AT g + (z)) (18) 

where AT+(z) is given by equation (fTKj) and z is the distance to the substrate (here AT+ > 
0, so that T — AT+ < T). Note the difference between the no slip boundary condition at 
the surface of the tip (if interactions between the tip and the monomers are weak) and the 
fact that the monomers on the substrate are attached. That latter condition is much more 
restrictive. A no slip boundary condition allows the monomers in contact with the tip to be 
renewed, which is not the case for the monomers in contact with the substrate if they are 
attached. For instance, if there was no specific interaction between the monomers and the 
substrate, one would have a no slip boundary condition also at the film-substrate interface. 
Then the slow aggregates encountered by the tip (Fig. 2.b) could rotate and reorganise 
within the film. According to our model, this difference results in the huge difference in 
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viscosity between the cases of suspended films and SI films, as illustrated by equations (fT2"j) 

and (flHjl. More precisely, the viscosity of the film can be calculated using equation (fTHj) 
and the bulk WLF law (or VFT law in the context of simple liquids) 

Vbuik(T) _ -Ci(T - T g ) 



50, 



511: 



log( 



Vbulk (T g 



T-T n 



(19) 



C\ and C*2 are two constants which depend on the considered polymer. The viscosity at T g , 
rj(T g ) is typically of order 10 12 Pa.s in the case of simple liquids. In the case of polymeric 
liquids, it is more appropriately defined a monomeric jump time of order 100 s |50j,|5l|. Close 
to T g , the slope of log(r](T)) is very large, such that an increase of T g of 20 K results in an 
increase of the viscosity up to 3 or 4 orders of magnitude depending of the liquids. When 
compared to the opposite case of a WI film, the same polymer films at the same temperature 
can have viscosities which differ by more than 6 or 7 orders of magnitude depending on their 
substrates. 



IV. ADSORBED POLYMER FILM 



We have just seen that when the film-substrate interaction is strong, the film is glassy 
when the size of the aggregates of slow subunits is comparable to the thickness of the film. 
This occurs when the correlation length of the 3D percolation problem is comparable to the 
thickness h of the film. The main point is that the monomers in contact with the substrate 
are "attached". As a consequence, a probe which moves within the film deforms slow 
aggregates, which thereby control the viscosity of the film. However, when the interactions 
with the substrate are finite, a monomer remains anchored on the substrate during a finite 
time. Thus the issue is: which interaction strength is required for the film to be glassy, or 
can we calculate the T g shift as a function of the interaction strength? Are typical adsorption 
interactions sufficient for inducing an increase of T g l Let us define e as the adsorption energy 
of one monomer with the substrate. If one could remove a monomer independently of the 
other ones, the life-time of monomer adsorption on one site would be: 

Td = T exp(^) (20) 

where To is a microscopic time of typical order of magnitude 10~ 12 s. For the film to be in the 
glassy state, the time r d must be at least equal to r g ~ 100 s. However, typical adsorption 
energies are of order a few T at most, which is far from required. 
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However, one must consider the number of monomers which must be removed collectively 
on the time scale r g . These are the ensemble of monomers which belong to the same aggregate 
of slow subunits and which are in contact with the substrate. One can see such an aggregate 
on Fig. 4. We call S = S(h,p) the number of such monomers. On a time scale smaller than 
T g , such an aggregate must be considered to be rigid. The time scale for removing such an 
aggregate from the substrate is therefore 

r, = r exp(^^) (21) 

where /(e) is the fraction of monomers in contact with the substrate which are adsorbed. In 
a first approach, this function can be considered to be of order unity and shall be omitted in 
the following. Of course equation (|21|) holds as long as the r.h.s is smaller than r g , since the 
aggregate does not persist on longer time scales. This equation breaks down at low e when 
the quantity eS/T is smaller than one. Then the time for removing the aggregates must be 
considered to be equal to r . Now, when considering equation (jUJ), one can understand 
why typical adsorption energies are sufficient for increasing the glass transition temperature 
of thin films. Indeed, when the size of the aggregates is comparable to the thickness h of 
the film, the number S is at least equal to iVf^ 3 , which is the number of monomers on one 
interface of a slow sub unit. Then, this number is of order a few tens, which can be sufficient 
for Td being comparable to r g , depending on the adsorption energy e. 

Let us describe now in more details the glass transition mechanism that we propose for 
thin films deposited on a substrate with interaction energy e per monomer. First let us note 
that the relevant temperature range is 

T-(h)<T<T g +(h) (22) 

where T~(h) and T+(h) are the glass transition temperatures of a freely suspended film and 
of a strongly adsorbed film and are given by equations (11) and (16) respectively. Indeed for 
T > T^ih) a film is fluid whatever the interactions with the substrate according to previous 
discussions. For T < T~(h) a film is glassy, even if the interactions with the substrate 
vanish, as in the case of suspended films. In terms of the fraction of slow aggregates, this 
temperature range corresponds to: 

p+(n) >p>p~{n) (23) 
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where p+(n) is given by equation © (n = h/£) and is larger than p 3 . , and Pc(n) is 
smaller than p 3D and is given by equation (|14p. Both curves p^(h) and p~^{h) are shown in 
Fig. 3. Both quantities correspond to a correlation length for the 3D percolation problem 
comparable to the thickness h of the film. One (Pc{h)) is below the percolation threshold 
(case of strongly adsorbed film), the other one (p^(h)) is above the percolation threshold 
(case of a suspended film). 

Let us start from the fluid state at T > T+(/i). By lowering the temperature below 
T^~(h) (i.e. by increasing p above p~(n)), the size of the aggregates increases and they 
attain a lateral extension Oi larger than the thickness h of the film: one enters the cross-over 
between 2D and 3D percolation. One can see such an aggregate on Fig. 4. We show in 
4fi| that this quasi-2D percolation problem with probability of occupation p and mesh size 
£ can be mapped, using a coarse-graining function $, into a 2D problem of mesh size h and 
probability of site occupation $^(p) with $^(p+(/i)) = p 2 c D . N is the number of iterations 
for transforming the film of thickness n on the scale £ (h = n£) into a 2D lattice of mesh 
size h. Using this procedure, one can demonstrate that the mass M' of an aggregate is [l^ : 

M'(h,p) « |$ (jV) (p+(n)) -$ (A °Go)r 72 « (|r 72M b+H -p|~ 72 (24) 

where the mass M' is given in unit of cubes of size h 3 in the renormalised system. On the 
other hand, the correlation length (expressed in the renormalised system, i.e. in units of h) 
is: 

Cj(M « \& N \pt(n)) - *W(p)|-« « (|)-^b?(n) -Pi"' 2 (25) 

where pjt(n) as a function of n is given in equation (JBJ). All these relations hold within a pre- 
factor of order unity which can be measured by numerical simulations The exponents 
2/2 ~ 1.333 and 72 ~ 2.389 are the 2D critical exponents for the correlation length and for 
the mass of the aggregates, respectively. In units of mesh size £, the mass reads 

M(h,p) w M\h,p)(jf-^ (26) 

where ^3 « —0.068 is the anomalous exponent of the 3D correlation function of the perco- 
lation problem. In physical unit length, the correlation length parallel to the film reads 

(\\~h(~)-^\pt(n)-p\-»> (27) 
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Note that the above result insures that for p = p~{n), i.e. when the 3D correlation length 
is equal to the thickness of the film, the size of the aggregates is equal to h as it should 
be. From equations (|24|) and (|26|) . one deduces that the number of monomers from one 
aggregate in contact with the substrate is (within a pre-factor of order unity) 

S{h,p) « iV 2 M(/ i ,p)(|)- 1 « N 2 (j) 1 ^(j)^ / *\pt(n)-p\-* (28) 

where N 2 ~ N% is the number of monomers on the interface of a subunit of size £. Equations 
((211), (J2EH2HI), together with equation (jHJ which gives the variation of p as a function of 
the temperature, and (jBJ) giving p+ (n), provide the variation of the mass of the aggregates, 
of the correlation length parallel to the film, and of the number of monomers from one 
aggregate in contact with the substrate, as a function of the temperature, for T < T+. At 
T = T~ all these quantities diverge. From equations (f2~Tj) and (J2EJ), one can calculate the 
time for removing an aggregate from the substrate. At the glass transition temperature of 
the film, this time must be equal to r g , which reads 

r d «roexp(^|^)=r g (29) 

This yields 

= (N 2 A^T\j) Mu3) H^\pt(n)-p\" 2 (30) 
J- ? ? T o 

This equation gives the fraction p of slow subunits in the film at the glass transition as a 
function of the interaction energy e. The glass transition temperature T of the film is then 
obtained using equation (jHJ) which relates the fraction p of slow subunits and the temperature: 
p — pl D = — j(T — Tg)/T g = —'yATg/Tg , where 7 = N^ 2 / ' F' \p z c D ) is positive. Using the 
relation p£(n) — p^P = —^/S.T g /T gi where AT~ < is the T g shift of thin suspended films, 
one obtains 

AT, _ 1 ^(f)'-- h ,„ A77 

T 9 - 7 ( rin(Vr ) ) ( «' + T, (31) 

or, using equation (JTUj): 

^^a-^k^;^.,) (32) 

T 9 7 f rin(r 9 /r ) 
This equation describes the cross-over between a suspended film and a film with strong 
interaction with its substrate, as a function of the adsorption energy per monomer e. At 
low values of e, one recovers the T g shift of suspended films, by setting e = in the above 
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equation. For strong interactions, the T g shift given by equation (}32|) is valid as long as it 
is smaller than AT 9 + , i.e. as long as the correlation length is larger than h. Equation (J3*2j) 
must be considered as leveling off at AT^ + when increasing e. Indeed, consider a system for 
which e can be tuned. At e = 0, the shift of T g is negative. The lateral size of the aggregates 
is infinite. When increasing e, the glass transition temperature increases. At e = e c , the shift 
cancels. At larger value, the shift of T g is positive. The lateral size of the slow aggregates 
has decreased but is still so large that the number of monomers of one aggregate in contact 
with the substrate is large enough for the desorption time of the aggregate to be equal to 
T g . In this regime the lateral size (\\ is still larger than the thickness h of the film. At some 
value e = e max , the glass transition temperature corresponds to a 3D size of the aggregates 
comparable to h. Then, increasing further the adsorption energy does not increase further 
the glass transition temperature since it would correspond to a correlation length much 
smaller than the thickness of the film. This plateau value is reached when e is a few e c as 
given below, and the corresponding T g shift is equal to AT+ . In any case our model predicts: 

AT-(h)<AT g <AT g +(h) (33) 



The interaction can be considered to be strong when 



> 1 (34) 



Tln(r 5 /r ) 

The interaction is weak, and the situation of the film is equivalent to that of a suspended 
film when 

( T1 7 , J 1/72 « 1 (35) 
T\n{T g /T ) 

Thereby, the cross-over between a freely suspended and a strongly interacting film of thick- 
ness h is given by 

Ec ~^ T iV 2 (|)^ (36) 

2/3 

With N c ~ 100, N2 ~ N c is of order a few tens. The time scale chosen for the glass 
transition is r g w 100s, which means that the quantity ln(r 3 /r ) ~ 30 typically. With 
h/l; = 10 for instance and [i ~ 0.5, one obtains then 

e c ~ 0.02T (37) 

This value for the interaction per monomer can be considered to be typical at the cross-over 
between negative shift of T g and positive shift of T g . Expressed in terms of the interaction 
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energy between the substrate and the polymer melt, one obtains e c ~ 0.4 mJ/m 2 . Note that 
the thicker is the film, the sharper is the transition between the two regimes (WI films and 
SI films), but the smaller (oc (h/£)~ 1 ' Vs ) the amplitude of the effect. One can write equation 
(E2J) as: 

TJh)-T-(h) 1 h * .. 

oc e 1/72 (38) 



T g (h) - T g (h) \ ( h l/v ^ lh2 . :j 



t 9 

to emphasize that (T g (h) — T~(h))/T g increases like e 1 / 72 ~ e - 418 i n the crossover region, a 
relation which may perhaps be tested experimentally. The glass transition temperature as 
a function of the interaction per monomer e is plotted on Fig. 5. Finally, the T g shift can 
be written as a function of the thickness of the film and of the interaction energy as: 

AT g (h,e) AT g (h) 



T T 

± 9 ± 9 



1 - I — 

e 



(39) 



in which AT~(h) is the negative shift in a non-adsorbed film (cf equation (fTTjO . The 

crossover value e c contains itself a h dependence and is given by equation (|3*B|l. Equation 

(JHSJ), together with equations (JTTJ) and (jSEJ), is the main result of this manuscript. 

As already discussed, equation (|3T)j) levels off at e = e max so that the T g shift is equal to 

AT^(h) in the regime of strong adsorption. One has thereby 

AT+(h) & 
ema X (h) = (1 - ^^)^c(h) = (1 + ^re c (h) (40) 

The ratio f3' / (3 is the ratio of the T g shift between the case of suspended films and that of 
strongly adsorbed films. Experimental data show that the typical shifts for 10 nm thick films 
are 20 K and 50 K in these respective cases. With 72 ~ 2.389 4f|, one obtains e max ~ 20e c , 
or given our prediction for e c for a film 10 nm thick, we obtain 

e max (h = 10nm) ~ 0.5T (41) 

which, expressed in adsorption energy per unit area, corresponds to about ~ 10m J/m 2 . 
Considering the respective values of e c and e max , one sees indeed that observing a negative 
shift of T g in a film deposited on a rigid substrate requires very well controlled surfaces, with 
very low interactions between the substrate and the monomers. 



V. DISCUSSION 



At this point, let us discuss some important aspects of the glass transition in thin films 
that can be deduced from our model. First, let us emphasise here that our assumption is 
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that the local dynamics on small scale is the same in the bulk and in thin films. The change 
of dynamical properties between both situations take place on larger scale due to percolation 
mechanism, which differ between the bulk and thin films. This assumption is underpined 
by the fact that the statistics are the same between both cases: t he p erturbations of the 



40j and do not extend 



thermodynamics due to the presence of interfaces is short ranged 
beyond 1 or 2 nm. 

An important consequence of our model is that the measured shift of T g should depend 
on the type of technique used for measuring the dynamics. This a consequence of the 
heterogeneous nature of the dynamics. The modification of the dynamics in thin films as 
respect to the bulk might be different for instance when considering viscoelastic properties 
and the diffusion of probes, as it has already been observed in the bulk, as a function of 



;emperature 



29 



30, 
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Note that the size of the probe plays a crucial role at this respect 



4l|. We expect also that the shift measured e.g. by NMR experiments, if any, should be 



different as that measured by mechanical experiments. Indeed, in a mechanical experiment 
as described in this manuscript, one probes only a part of the dynamical spectrum (as defined 
by equation (jHJ)), according to the model proposed in j^Q, whereas a NMR experiments 
probes a more local environment. 

However, the T g shift should also depend, for a given technique, of some precise aspects of 
the experiment. For instance, we argued that a suspended film is more fluid when considering 
pulling a probe in the film, while for the same experiment we expect a strongly adsorbed film 
to be more viscous. Consider now an experiment in which one shears the film between two 
rigid substrates, which have no interaction with the considered liquid. Then, the T g shift 
should be different from that measured by pulling the probe and discussed in this paper. 
Indeed even at temperatures above T~ the slow aggregates are strongly anisotrope and are 
very elongated in the direction parallel to the film, if T is not too far above T~ . Then, in a 
shearing experiment as discussed here, one needs to deform them given the no-slip boundary 
condition. Some dilatancy is also expected, depending on the shearing rate. In other words, 
we expect the viscosity of the film to be strongly anisotropic. 

Another issue is a more precise calculation of AT+. Indeed, our result here (equation 
(115(0 has been obtained only within a pre-factor, assuming the relation p1 D —p~ ~ {h/^)~ 1 ' u , 
so that a fraction of order one of slow subunits on one interface are connected to the other 
interface by a slow aggregate. But we did not specify the value of this fraction. Thereby, 



17 



the pre-factor (3' in equation (JT5|l must be considered to be an adjustable parameter, which 
we choose to be consistent with a T g shift of about +50K for str ongl y interacting films 



of thickness lOnm, which is close to experimental results [ill, |13|, Il4l ]. For the type of 
experiments discussed here (e.g. probe motion, Fig. 1 and Fig. 2), the viscosity of the 
suspended film is controlled by the slowest percolating subunits, leading to an equation 
analogous to equation © with p^(n) instead of p z c D . The percolation threshold in the 
direction parallel to the film is a well defined state, which is also the case in the bulk. It allows 
thereby to calculate the shift of T g between suspended films and the bulk as corresponding 
to the shift of the percolation threshold. In particular we expect the glass transition in thin 
suspended film as discussed here to be similar to that in the bulk. The situation of strongly 
adsorbed film is very different. The fraction of slow subunits at the free interface connected 
to the lower interface is never zero, but decays exponentially with the thickness h of the film: 
the percolation between both interfaces is not a well defined state, but occurs continuously. 
In addition to the broad distribution of relaxation times already present in the bulk and in 
thin films, this fact tends to increase the width of the glass transition, as compared to the 
bulk or as compared to suspended films. This fact might be observed experimentally. For 
our purpose here, it makes a calculation of the T g shift more involved, as it requires a more 
precise model as that presented here for calculating the viscosity of the film as a function of 
the statistics of the slow aggregates. 



VI. COMPARISON WITH EXPERIMENTAL RESULTS 

We have described in this paper a mechanism which leads to a T g shift in thin polymer 
films. This shift depends on the film thickness and it can be either positive or negative 
depending on the interaction between the monomers and the substrate. Two limiting cases 
must be considered when dealing with equation 1) the case of very weak interactions. 

Then, the T g shift is negative and equal to AT g {h). The variation of AT g as a function 
of h is then given by a power law AT g oc h- x l v \ 2) the case of strong interactions. When 
increasing e, equation (jHUj) must be considered as leveling off when e becomes larger than 
tmax which is a few tens of e c , so that the T g shift saturates at This is a direct 

consequence of the model proposed here. Increasing indefinitely the interaction energy does 
not change the T g shift. In the strong interaction limit, the T g shift is equal to AT+(h) and 
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has the same power law dependence on h: AT g oc h x l Vz . In the crossover region, the T g 



shift is a function of both h and e. Since e c itself depends on h, the T g shift is not exactly 
a power law of h in this regime and the complete equation ([39)1 must be used to fit the 
data. When considering experimental results, one needs thereby to check whether one deals 
with one of the two limiting case, or whether the system is in the cross-over in between. Let 
us now discuss some experimental data regarding the shift of T g on various substrates and 
conditions within this framework. 

Kawana and Jones Q have reported various measurements of the glass transition temper- 
ature T g of polystyrene (PS) films deposited on silicon wafers by spin coating, as a function 
of the thickness h of the film. T g is measured from thermal expansivity curves obtained by 
ellipsometry. In these experiments, T g is shifted downwards, which corresponds to a case 
where the interactions between the polymer film and the substrate are weak. The data were 
fitted using the empiric relation 



with a best fit value 5 = 1.1. This value is indeed very close to our prediction (eq(ll)). 
However, the dispersion of the data is too large for discriminating precisely the value of this 
exponent. The amplitude of the observed T g shifts corresponds to a length A = 0.8 nm. We 
shall compare this value to our prediction for a film on a non-adsorbing substrate: 



with /3 = -fxF'(pl D )/Nc /2 ~ 1/3,/ . The values /i « 0.5, F'(p 3 c D « 0.11) « -4, N c « 100, 



v ~ 0.88, a ~ 0.6 nm give A = a(3 v ~ 0.7 nm, which is in remarkable agreement with 
the measured value. Our prediction is thereby in quantitative agreement with Kawana and 
Jones data, and we conclude that their experiments correspond indeed to the case of a non- 
interacting film, for which the effect of the interactions with the substrate is negligible. We 
retain here that the typical T g shift for a non-interacting film of thickness h = 10 nm is 
about 20 K (for a bulk T a of the order 375 K), or about 10 K at h = 20 nm. 

n 

Hartmann et al |Xf| have reported dielectric measurements in isotactic poly(methyl 
methacrylate) (i-PMMA) deposited on aluminium electrodes evaporated on glass substrates, 
which show a decrease of T g as h decreases. The shift is of the order 10 K for h = 20 nm, 
which in this case also coincides with the value mentioned above for a negligibly interacting 




(42) 




(43) 
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film (with here T g of the order 330 K). We deduce thereby that in that case the interactions 
were indeed negligible. In another series of experiments, they have also performed ellip- 
sometric measurements in i-PMMA deposited on silica wafers, which show an increase of 
T g as h decreases. This apparent contradiction is interpreted as due to a difference in the 
polymer-substrate interaction strength, because the experiments have been performed on 
different substrates. Hartmann et al interpreted these measurements in terms of a model by 
Herminghaus 2, which gives an expression for the T g shifts as a function of the ratio a/E 
(a the surface tension, E the Young's modulus of the polymer at T g ). However, within Her- 
minghaus model, only decreases of T g are predicted, disregarding of the polymer-substrate 
interaction. Changing the sign of the T g shift as predicted by Herminghaus for interpreting 
' iartmann and co-workers experiments is thereby an unsubstantiated extension of this model 
. On Fig. IHlwe plotted Hartman and al's results obtained by dielectric measurements, for 
two different molecular weights. The straight line represents the extrapolation of the data 
according to equation (}43|) . with a slope — I/V3 = —1.136. The data are compatible with 
this slope, even though they are also not precise enough to test this exponent. The ordinate 
at log h = gives the coefficient A m 1 nm. Equation (f3H|) of the present paper provides 
an explicit dependence of the T g shift as a function of the interaction strength e in which 
AT^(h) is the negative shift obtained for a non-interacting film. As discussed above, the 
crossover value e c contains itself a h dependence (cf equation (}3*HjO and in this regime, the 
T g shift is not precisely a power law. Nevertheless, if the /i-dependence in e c is neglected in 
equation ()39j) . the intersection of the solid line in Fig. with the ordinate axis provides the 
pre-factor of equation (jSHJ), i-e. log((e/ 'ec) 1 ^ 2 — 1) ~ 0.7, or (e/e,,) 1 / 72 « 6, that is e « 75e c 
-assuming log A = 0, i.e. A fa 1 nm. From that we deduce that these experiments have 
been performed in the regime of strong interactions. In Fig. d) we have indicated the slope 
—1.13 which fits reasonably well the data. Fryer and co-workers Q| have measured the 
T g shift of PMMA films of various thicknesses deposited on silicon oxide, using a thermal 
probe technique. Here also, positive shifts have been reported, which correspond thereby 
to non-negligible interactions between the polymer film and the substrate. However, the 
positive shift here is not as large as in the previously discussed case. In Fig. |HJ), the ordinate 
at log/i = gives log((e/ tc) 1 ^' 2 — 1) ~ —0.2, that is e fa 3e c (with A fa 1 nm). We are 
thus here in an intermediate regime, where the ratio e/e c is not very large with respect to 
one. Therefore the variation of AT g versus h should be better described using the exact 
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TABLE I: Values of the T g shift for a thickness h = 20 nm, for PS and PMMA films deposited on 
substrates with various interfacial energies, from the data in Fryer et al (2001). For PS: T g = 375 
K; for PMMA: T g = 392 K. The values given for e/e c are based on a shift of 10 K for h = 20 nm, 



and the h variation in e c is neglected. 

Polymer a (mJ/m 2 ) AT g e/e c 

PS 0.8 -8 0.02 

PS 1.7 -4 0.3 

PS 6.5 +30 27 

PMMA 0.5 -8 0.02 

PMMA 1.3 -5 0.2 

PMMA 3.3 +11 6 



equation However, in the above discussion, the data have been extrapolated using the 

power law h~ 1A36 . Indeed, the data are not precise enough and the h range is too small to 
discriminate the strong adsorption limit and the cross-over regime, as far as the variation 
of T g versus h is concerned. The amplitude of the shift, however, allows for this discrimi- 
nation. In another series of experiments, Fryer et al 17] have measured the T g shift in PS 
and PMMA films deposited on substrates with tuned interfacial energy. This was achieved 
by using as substrates octadecyltrichlorosilane (OTS) films deposited on silicon wafers and 
exposed to various doses of X-Rays. The authors claim that, at a given thickness, the T g 
shift varies linearly with the interfacial energy a, which at first sight may seem in apparent 
contradiction with our prediction (equation (|39|) ). However, we stress that the adsorption 
energy e and the interfacial energy a are two different quantities, and that the relation 
between these two parameters is not straightforward \v\ . . Nevertheless, the data pub- 
lished by Fryer and co-workers offer the opportunity to estimate the ratio e/e c using our 
model in the different cases considered in [17J. For each polymer (PS or PMMA), a series 



of results obtained for 3 different thicknesses are presented, for 3 values of the interfacial 
energy (obtained by an X-Ray irradiation with 3 different doses). The thickness-dependent 
T g shift is either negative or positive, according to the interfacial energy. The value a c at the 
crossover is obtained as the point where AT g = on each curve AT g versus a, for each given 
thickness h. It is observed that all curves (obtained for various h values) intersect at a point 
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of coordinates AT g = 0, a ~ 2mJ/m 2 . In terms of interfacial energy a, the crossover should 
therefore correspond roughly to the value a c ~ 2mJ/m 2 . Note however that the surface 
tension a is not the relevant parameter for describing the strength of the interaction as far 
as the T g shift is concerned: the relevant paremeter is the adsorption energy per monomer. 
More precisely the relation between the T g shift and a depends on the considered substrate 
and polymeric liquid. We expect only the relation between DeltaT g and the adsorption 
energy e to be universal as given by equation (39). Using equation (j3Uj) . and by taking 
/ST g (h = 20nm) — 10K, one can deduce the ratio e/e c . The results are summarised in 
Table I. One can see there that the regime of strong interaction is reached when e/e c ~ 27. 
Note that the T g shift of 30 K for a film 20nm thick is in the plateau regime. Thereby, our 
estimate here of e/e c is an underestimate, and the real adsorption energy might be higher. 



VII. CONCLUSION 



In conclusion, we have proposed here a model for calculating the T g shift of thin polymer 
films, as a function of the adsorption energy between the polymer and the substrate. We 
show also how to calculate the corresponding viscosity of the film using the bulk WLF 
law with the corresponding temperature shift. Our prediction could be tested directly by 
pulling a tip immerged in the film, such as an AFM tip, or a probe pulled using an optical 
tweezer for instance. As a consequence of this model, we expect the viscosity of a suspended 
film to be strongly anisotropic. For instance, the viscosity measured by shearing the film 
between two susbtrates (without interaction) should be much larger than that measured 
by pulling the probe as discussed in this paper. Such experiments could be performed 
in the two limiting cases of non-interacting films and strongly interacting films, and also 
in the cross-over regime of intermediate adsorption energy. However, testing our general 
relation quantitatively would require a precise knowledge of the interaction between the 
monomers and the substrate. Finally, our prediction of long ranged changes of viscoelastic 
properties should be of importance for explaining the reinforcement of filled materials such 
as filled elastomers 13, 1^ and should help in tailoring their properties by adjusting the 
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adsorption energy between the particles and the matrix. 
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VIII. FIGURES 
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FIG. 1: We represent here a body floating in a thin suspended film, seen from above. The body 
can be thought of emerging from the film on both interfaces. We consider in Fig. l.a a film at 
the bulk T g . Since the film is suspended, the slow subunits in Fig. l.a, which do not percolate in 
the direction parallel to the film, can move freely in the film to open the way to the probe: the 
viscosity is smaller than in the bulk. By lowering the temperature (b), the slow subunits percolate 
in the direction parallel to the film. They now control the viscosity of the film as measured by the 
motion of the probe. 
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FIG. 2: We consider here a film at a temperature T higher than the bulk T g . In Fig.2.a, the 
tip is at a distance larger than the size ( of the slow aggregates. The slow aggregates can move 
freely around it, and the viscosity is smaller than the viscosity at T g . In Fig.2.b, the tip is closer 
to the substrate. Moving parallel to the film requires deforming the slow aggregates: the viscosity 
is larger than that at the bulk T g . 
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FIG. 3: p£(h) and Pc(h) as a function of /i. Equivalently, this graph provides the 3D correlation 
length of the percolation problem -within a prefactor-, as a function of p. A freely suspended film 
of thickness h is glassy at temperature T if p (related to T by by equation (8)) is larger than p+(h). 
A film with strong interactions with the substrate is glassy if p is larger than Pc{h). Thereby, at 
any temperature lower than T g + AT~ a thin film is in the glassy state, while at any temperature 
higher than T g + AT^ 1 " the film is in the liquid state. The cross-over between low interactions and 
strong interactions corresponds to the region Pc(h) <p< p£(h). 
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FIG. 4: An aggregate of slow subunits in a thin film, in the cross-over temperature range i.e. 
T g + &T g < T < T g + AT 9 + . The lateral extension of this aggregate is Cy, and the number of 
subunits it contains is M(h,p) 



29 




FIG. 5: Variation of the T g shift of a thin film, as a function of the interaction energy between 
a monomer and the substrate. The crossover between suspended films and films with strong 
interactions with the substrate takes place at e = e c ~ T. 
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FIG. 6: Data from dielectric measurements in i-PMMA films deposited on aluminium by Hartmann 
et al (2002). T g decreases as h decreases, so that AT g = T g (h) — T g is negative. Two different 
molecular weights are shown (□ : M w = 44900 g mol -1 , T g = 330.6 K, ■ : M w = 164700 g mol -1 , 
T g = 331.2 K). The line has a slope -1.136. 
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FIG. 7: Data from ellipsometric measurements in i-PMMA films deposited on silica wafers by 
Hartmann et al (2002). T g increases as h decreases, so that AT g is positive here. This case 
corresponds to strong interactions between the polymer film and the substrate. (M w = 44900 g 
mol -1 ). The solid line has a slope -1.136. 
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FIG. 8: Data from thermal probe measurements in PMMA films deposited on silicon oxide by 
Fryer et al (2000). T g increases as h decreases, so that AT g is positive here. This case corresponds 
to intermediate adsorption energy between the polymer film and the substrate. (T g = 392K). The 
solid line has a slope -1.136. 
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